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We review several results that have been obtained using lattice QCD with the staggered 
quark formulation. Our focus is on the quantities that have been calculated numerically 
with low statistical errors and have been extrapolated to the physical quark mass limit 
and continuum limit using staggered chiral perturbation theory. We limit our discussion 
to a brief introduction to staggered quarks, and applications of staggered chiral pertur- 
bation theory to the pion mass, decay constant, and heavy-light meson decay constants. 
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1. Introduction 

Currently, a large number of lattice res ults have appeared using staggered fermions 
to discretize the quark fields.^ ^ ^ '^ l ^ l ^ l ^ l What contrasts these calculations with ear- 
lier lattice calculations is the use of 2-1-1 dynamical flavors of light quarks, with up 
and down quark masses down to ms/10 (where the lattice strange quark mass is 
roughly its correct value), multiple lattice spacings, and the ability to understand 
the light quark and continuum limit. This last step allows one to minimize system- 
atic errors in many of these results, so as to make an accurate extrapolation to the 
physical values of the light quark masses. 

For staggered fermions, this step is highly non-trivial. Introductions of scaling 
violations which are rather large at finite lattice spacings cause strong effects. These 
violations are reduced using improved actions, but are still not negligible .[^Because 
many of these violations add to the squared pion mass, it is perhaps not surprising 
that when attempting to construct the relevant chiral perturbation theory, they 
must be included as chiral symmetry breaking terms, similar to the way the quark 
masses enter. 

This review is intended to be focussed on this aspect of analyzing lattice data. 
We begin by outlining lattice QCD with staggered quarks in Sec. [21 We continue 
by discussing the general procedures for deriving a low-energy effective theory for 
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QCD with staggered quarks in Sec. [3l We then use this to calculate some pion- 
related quantities in Sec. [4l and include heavy quarks in Sec. [5l Overall, we will 
focus primarily on results from chiral perturbation theory that have been applied 
to lattice data to obtain relatively small systematic errors as a result. There are 
many other quantities that have been studied, and we will briefly touch upon them 
in Sec.m 

2. Lattice QCD v^rith Staggered Quarks 

We begin by formulating Quantum Chromodynamics on a discrete spacetime lattice 
with a uniform lattice spacing a. The quark fields live on sites, and the gluon fields 
live on links, although we will not discuss the gluons in detail. While focussing on 
the quarks, we will discuss everything in terms of the free theory for simplicity, but 
the results carry over to the interacting case. 

The simplest discretization is given by replacing the derivatives in the free Eu- 
clidean quark action by difference operators: 



'5'quarks — ^ ^ ^ 

f.x L 



(1) 



^ ^ lfiQf,x {Qf,x+afi — Qf,x-afi) + mfQf^xQf,x 

where / runs over the three lightest quarks u, d, and s0 a; is a spacetime point, /i is a 
unit vector in the ^-direction, and to/ is the mass of the quark flavor /. This action 
has the well-known doubling problem, which can be seen in the momentum-space 
propagator for a given flavor / 



~«E^7Msin(ap^) -ham/ 



V • (2) 

This has a low-energy mode as a — > at p = 0, but also when any component of 
p is near 7r/a or — at one of the corners of the Brillouin Zone. The app eara nce of 
these additional 15 doublers is a general result of discr etizing fermion s.l^^^ There 
several techniques to rid oneself of the doubling problem j-*^*-* l ^^ l '*^^ l '*^3 | 14 | g^j^^j ^-t^w 
focus here on the staggered discretization.HH 

To arrive at the staggered quark formalism, we realize that when going to a 
discretized Euclidean space, one breaks the 5*0(4) rotational symmetry down to 
a discrete subgroup composed of finite rotations by 7r/2 (denoted by SWi). This 
implies that we can perform a redefinition of the quark fields in Eq. ^ that di- 
agonalizes the 7 matrices, thereby decoupling the different spinor components of 
Qf. We then keep only one spinor component per field and we have decreased the 
number of degrees of freedom by a factor of four, and giving us the action 



r — «4 \ " 

■^quarks — ^ ^ ^ 



2a 



(3) 



^The heavy quarks play no dynamical role and can be integrated out at the scales usually used in 
lattice simulations. 



I 



I 
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where rj^^x = X]i/</i(~l)^'^ > ^^'^ ''Pf single-component field remaining after the 

field redefinition and dropping the other spinor components. To show that this set 
of 16 single-component fermions (per flavor /) corresponds to four four-component 
Dirac fermions is a straightforward exercise that we will not pursue here. The most 
important point here is that different species (or "tastes") of a given flavor are 
identified (in momentum space) with different corners of the Brillouin Zone."'^'' 

It is possible to take only a single flavor and add a more general mass matrix 
which can split the degeneracy in Eq. ([3]) so as to describe four-flavor QCD (the 
four tastes are identifled as four different flavors). However, this procedure has 
technical difficulties making it rather tough to implement numerical. What is done in 
practice is to create one staggered field for each flavor one wishes to simulate. After 
evaluating the path integral over the fermions, we obtain the fermion determinant 
for each staggered fleld, which is a determinant describing four tastes. Since these 
four tastes are degenerate in mass, the fourth-root of this determinant should then 
be a determinant describing a single taste of the staggered quark.^I^ 

This "rooting" technique is rather common in staggered simulations, but it is 
not without theoretical issues. At finite lattice spacing, there are interactions among 
staggered quarks which violate the taste symmetry: the symmetry which rotates the 
tastes among themselves and is a symmetry of the free theory. In the interacting 
theory, the exchange of high momentum gluons can change the taste of a quark. 
This is shown in Fig. [1] and one can see that by exchanging a momentum of order of 
the cutoff (g « tt/o in the figure) does not move the quark far off shell, as it would 
with other discretizations, but instead changes the quark to a different taste with 
a small amount of momentum. 




Fig. 1. A four-quark diagram which can introduce taste- violating interactions due to the exchange 
of high-momentum gluons. q here is assumed to be of order it /a. u and u' are up quarks of different 
tastes, as are d and d' . 

The question remains as to whether or not it is valid to take the fourth root of the 
staggered determinant before one takes the continuum limit. If these two operations 
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did not commute, then a rooted-staggered theory would not be equivalent to a 
one-flavor continuum theory in the continuum limit. Although there is no proof 
that exists that proves the validity of this procedure, there are indications that the 
rooted theory is legitimate. There is a large body of supporting evidence in favor 
of the v alidity of the ro oted theory, as well as some evidence implying this trick is 
not valicP^^l^^I^^^QllII ('^g refer the reader to Ref. [22] for an excellent review of the 
current status of the rooting procedure). Assuming the validity of this trick, we will 
now pursue the application of XPT to the rooted staggered theory. 

Although the rooting procedure gives the correct number of degrees of freedom in 
the sea sector, this does not mean that we cannot see the effects of these additional 
tastes. In fact, we will see that staggered chiral perturbation theory, even after 
correcting for the fourth-root trick, shows a significant dependence on the different 
tastes of quarks that arise. 

3. Staggered Chiral Perturbation Theory 

Since one cannot simulate at physical u and d quark masses currently, one must 
extrapolate to the physical masses from the larger masses actually used. For this 
one needs chiral perturbation theory (XPT). However, continuum XPT is not ap- 
propriate for a lattice theory, as there are other sources of chiral symmetry breaking 
besides the mass term. For staggered quarks, this additional chiral symmetry break- 
ing arises primarily from the taste-violating interactions discussed in the previous 
sectiorlB. As such, one must appropriately modify the c hiral theo ry to take into 
account these additional chiral symmetry breaking effects, ^^^^^3125] construct a 
staggered chiral perturbation theory (SXPT). 

The first step in this procedure is to write down the Symanzik effective actioiJ^ 
for the lattice theory. The idea is that, given a lattice action S'lat, one can write it 
as an expansion in a about the continuum limit 



where Co is the continuum QCD Lagrangian. For the staggered quark and gluon 
actions, this expansion is in a? (there are no errors linear in a), so only even values 
of n are included in this sum. We will work only to 0{a?) in this review. 

The first term, Cq^ is straightforward, as it is merely the continuum Lagrangian 
corresponding to the lattice theory. For higher-order terms, all possible operators 
that satisfy the lattice symmetries must be included, making this procedure pro- 
gressively more difficult when working beyond £2 • We will discuss the Lq and some 
of the £2 terms in detail, although the £4 terms have been worked out as well.^^Sl 
In current lattice simulations with staggered quarks, we have a^Ag^^^ ~ mJ/A^ 
with A^ the chiral scale. Since in the continuum chiral expansion, ~ mj, we 
use a power counting scheme such that ^ a? as well. This implies that we will 

''There could be other forms of discretization errors but taste violations are dominant. 
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omit terms in £2 that have powers of the quark masses or derivatives, as these 
terms would necessarily be of order 0{p^a?) (we use p generically for any power of 
derivative or pion mass for simplicity), which is one higher order than that which 
we are working to. 

Once we have the Symanzik action worked out to a given order, we can then 
transcribe the Symanzik terms to t he und erlying chiral theory. This process is well 
understood for the continuum theor j^^ZEHl ^nd we will work through the basic details 
when discussing Cq, as going through this process for the higher-order terms is a 
straightforward extension. 

3.1. Co 

Although the Cq term is simply the continuum Lagrangian, recall that this action 
is one with four tastes for every flavor of quark. For QCD with N flavors, we have 

Co = qL{x)pqL{x) + qR{x)pqR{x) + qL{x)MqR{x)+qji{x)M\L{x) + • • • , (5) 

where the • • • refers to terms which do not include the quark fields, such as the pure 
gluon Lagrangian. The quark fields are 4Af-component objects with the 4 tastes 
for each of the N flavors, and M is the AN x 4iV mass matrix (note that the four 
tastes for each flavor have the same mass). Also, the mass matrix is real, but the 
forthcoming spurion analysis is simpler by considering it to be a complex operator; 
in the end we set M = . Cq has an SU{4:N)l x SU('iN)ii chiral symmetry in 
the limit M 0, where, (with L, R e SU{4:N)l,r) 

qL LqL, qR RqR (6) 

is a symmetry of Eq. Q . The chiral symmetry is broken down to the vector sub- 
group SU{A:N)y by a non-vanishing quark condensate, which gives rise to (4A'^)^ — 1 
Goldstone bosons. These bosons are encapsulated in a field E = exp [i<I>//], where 

E ^ LEi?t ^ (7) 

and $ is the 4A^ x 4A^ matrix which contains the Goldstone bosons. For the case 
where A'^ = 3, we have 

/ U TT+ K+\ 

<P= in- D , (8) 

Va'" S J 

where U = Y^lti UtTt, etc., and Tt € {^5, «^m5i ^Cm'^i ^^}- We use the Euclidean 
gamma matrices with = ^^^1/, ^^5 = ^^^5, and ^/ = / is the 4x4 identity 
matrix. We write Eq. ([8]) in the flavor basis, where U ~ uu, D ~ dd, etc. 

Since the quark masses are non-zero but small compared with Aqcd, we can 
include them as a soft breaking of the chiral symmetry. To do this systematically, 
we promote the mass matrix to a spurion field with the transformation 

M LMR\ RM^L'' , (9) 
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which makes Cq invariant under the chiral symmetry group. Thus, using S, S^, M, 
and , we can write down a low-energy effective action for Cq (after setting M = 



This is just a continuum XPT expression for AN flavors of quarks in Euchdean 
space. / and fi are unknown parameters that are directly related to the tree-level 
pion decay constant and chiral condensate. 

The methodology of determining the low-energy chiral Lagrangian from the 
Symanzik action is the same for C2, ^4, and so forth. We need to know the operators 
that exist at the level of the Symanzik action, which can, in general, break the chiral 
symmetry of the massless QCD Lagrangian. By promoting objects that appear in 
the Symanzik action to spurion fields with specific transformation properties (as we 
did for the mass matrix), we can make the theory invariant and then from this, we 
can generate the corresponding operators at the chiral level. 



Terms in the Symanzik action contributing to £2 come from several sources. They 
can only arise from two-quark operators or four-quark operators, but we can im- 
mediately drop the two-quark operators from our analysis here. This is because a 
two-quark operator has dimension three, which goes to dimension one when we add 
in the required factor of a^. This means we need three powers of mass, in the form 
of the quark mass or derivatives, and then the overall term will be at least of order 
C(a^p^), which is higher order than that to which we are working. The four-quark 
operators are the only terms we need currently (one can easily see that higher num- 
bers of quark fields will only contribute at higher powers of a^). Effectively, we are 
looking at the terms in £2 in the limit — and M — 0. 

In determining the operators that arise in C2, one must first enumerate all of 
the operators that are consistent with the full lattice symmetry group at finite a. 
The full symmetry group was discussed for = 1 in Ref. [23] and for general N in 
Ref. I25l For non-zero quark mass (assuming n of the quarks are degenerate and m 
are non-degenerate, with N = m+n), this symmetry group is ?7(n)vec x (C^(l)vec)'" x 
r4 X S'W4,diag- Here, C/(m)vcc and [/(l)voc are flavor number symmetries, r4 is the 
Clifford group with generators and S'W4,diag C 5*0(4) is the group of hypercubic 
rotations in Euclidean space. 

The number of operators that contribute to £2 is around 25, where spin and 
taste indices are uncoupled (there are 10 operators where the spin and taste indices 
are coupled, but these do not contribute at this order). Many of these, however, do 
not need to be analyzed, as it was shown in Ref. [25] that only "odd-odd" four-quark 
operator^ need to be included in our list. With this realization, we can take the 

'^Odd bilinears are those where the staggered quark fields are separated by an odd number of Unks 



Mt): 




(10) 



3.2. £2 
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analysis for the N = 1 case from Ref. [23] and make the replacement of the taste 
matrices 

6 - ^ = 6 ® Iat (11) 

where Iat is the N x N identity matrix in flavor space, and t can be any of the 16 
tastes. 

Although we will not go into detail on the analysis for all of the operators, 
we will sketch the steps for a single set of operators. For simplicity we will pick a 
four-quark operator that has a vector spin structure and arbitrary taste structure: 

= ^ [qn (tm ® ^) qR + U (it. ® T^f ^) qL] ' , (12) 

where the g's and g's are the same 4iV-component objects above, and the taste 
matrices are going to be set equal in the end, T^'^ = T^^^ — T^^\ r(^) can either 
be a pseudoscalar (^5^') or tensor taste, since this has to be composed of 

odd bilinears.^23l We give these taste matrices the spurion transformation under our 
SU{AN)l X SU{'iN)R chiral symmetry: 

Tr^RTrR\ Tl^LTlL^ (13) 

to make Eq. (|12p a chiral invariant, just as we did for the mass term in Cq. Since 
the only objects that appear in this term are the taste matrices and the quarks, 
we construct the chiral-level operator from Tr l, S and S^. The only operator that 
arises is 

^ Tr[TLETflSt] . (14) 

This term leads to two different operators when we set T^^"* = T^^^ = T^^\ with 
T^^) either the pseudoscalar or tensor tastes: 

O- ^ TrKf ^Sef 'st] , 

tj,<i' 

Working through the rest of the operators in Ref. [531 we get a total of eight 
operators at the chiral level. We write these as 

a^Vx = a^U + a^W (15) 



on the lattice. 
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^U^C, Tr(ef )s^f )Et) + C,\ 5:[Tr(e^Sff )£) + h.c] 

= C2v\ E[Tr(er^S) Tr(e^S) + h.c] + C,a\ E[Tr(C^f S) TrC^^f ^E) + h.c] 

+C5V^E[Tr(^r^S)Tr(^Wl]t)] +C5Ai^[Tr(^(f S)Tr(4::)l]t)] . (17) 

Thus, the full Lagrangian at the chiral level including both £o and Li in Eu- 
clidean space is 

C<- = i- Tr (9^Ea^Et) - ^ Tr (AfE + Y) M) + — ^ Tr[$]2 + c?V , (18) 
8 4 3 

where we have dropped terms that come from higher order in the joint m^,a^ 
expansion. We have also included here the term that is required by the anomaly 
(note that Tr[<i>] is merely a sum of the flavor- neutral taste-singlet mesons), and 
we refer the reader to Ref. [28] for more details. Gasser & Leutwyler have worked 
out the terms that come in at 0{p'^, m^p^, m''j^and Sharpe & Van de Water have 
included the terms which arise at 0{p^a^, m^a^, a'* ),|29lbut we will not include these 
contributions to the action explicitly here. They will appear as analytic terms when 
we calculate quantities to one-loop order, however. 

4. Pions 

We can now use Eq. (fT8|) to calculate light meson (referred to from now on gener- 
ically as "pions" ) quantities. We first expand the field E to quadratic order in the 
$ field, and we can determine the tree-level masses for the pions. We get 

mly^t = + my) + a^A(^t) (19) 

with t the taste index, and x and y are any of the quarks. For a = we get the 
standard relationship between the pion mass squared and the quark masses. The 
additional term depends on the taste, t, and we have 

: 

: i| (Ci 3C3 + C4 + 3Ce) 
f 

: i|(2C3-f 2C4 + 4C6) 

^ (4C3 + 4C4) . (20) 



A(e5) 


= Ap 


A(^;.5) 


= A^ 


A(^;..) 


= At 


A(?,) 


= Ay 


A(6) 


= A, 
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Note that, to 0{a^), there is a remnant SO (A) symmetry that keeps certain tastes 
degenerate, and so we use labels that refer to the representations under S'0(4) 
rotations. This was first seen for the iV = 1 case in Ref. [23|and carries over to the 
general TV casePS Also, only coefficients from U contribute to the tree-level masses. 
W does contribute here as well, but by giving additional two-point vertices much 
like the anomaly does. It only affects the flavor-neutral (diagonal) pions which have 
vector (V) or axial-vector (A) taste. Expanded to quadratic order in the fields, we 
have 

a^U'^^^U^ + D,, + S^ + ---f + V^A (21) 

where 

S'y = ^{C2v - C5v) , (22) 

and similarly with A. These are of the same form as the anomaly term, and so we 
treat them in the same manner. The W terms for the V and A tastes as well as the 
anomaly term for the singlet tastes all act in such a way to add off-diagonal terms 
to the mass matrix. By diagonalizing this matrix, we can write everything in terms 
of physical fields, 7r° , 774 , 77^ |j Of course, in the end, we send toq 00 since the 77^ is 
heavy on the chiral scale, but we cannot do this for the ry^ and rj'y. 

The procedure for diagonalizing the mass matrix, or equivalently calculating the 
full propagator, for the flavor-neutral flelds is described quite well in Ref. [30l and 
we just state the results here. For comparison, we first write the propagators for the 
flavor-charged fields with taste t, or fiavor-neutrals with taste T or P: 

1 

r + m-y^t 

For the fiavor-neutrals, the propagator can be a little more complicated. For the 
taste- vector, taste-axial, or taste-singlet we have 

with for the axial and vector tastes: 



G.y.til') = 2^ 2 ■ (23) 

_L m ^ 



G^yM)- .^. (24) 



('?^ + ™L,v(A))('?^ + ™Lv(^)) 



(g' + ^L,y(A))(g' + ^L,v(A))(g' + "^L,v(A)) 

(g' + "i^o V, .0(g2 + vMi)(g' + K 



(25) 

We have written Eq. ([^5]) explicitly for = 3, but one can see the pattern. There 
are N factors of (q^ -I- m^) in both the numerator and denominator. Those in the 

"^We refer to these as physical because for the singlet case (relevant in the continuum) this amounts 
to changing from the flavor basis: U,D, and S, to the physical basis: n^,ri, and 77'. Of course, the 
axial and vector taste fields are not strictly "physical," but we will use the same terminology since 
they all result from diagonalizing the flavor-neutral mass matrix. 
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numerator correspond to masses of the pions in the flavor basis and those in the 
denominator are those in the physical basis. The result is the same for the taste 
singlets, but we can take toq — > oo in those terms to get 

One final issue we have avoided thus far is the fact that we have an incorrect 
number of tastes per flavor in our theory. Our theory has AN species of quarks, or 
four tastes for the N flavor s. The simplest procedure, for one-loop calculations, is 
to use a quark-flow analysiJ^^^ to determine which quarks are "sea" quarks and 
thus correspond to quark loops. These are the loops which are removed by taking 
the fourth root in lattice simulations, and we must correct for this in the chiral 
theory. Since each quark loop can have four tastes, we merely multiply these terms 
by a factor of 1/4 to reduce the theory to the correct number of tastes per flavor. 

The multiplication factor arises in two places. Terms which are connected (only 
include the -|-m^) parts of the propagators) and those which are disconnected, 
coming from the V portion of the flavor-neutral propagators. For the latter case, 
the correction is rather straightforward. As explained in detail in Ref.1251 we change 
the values of the parameters 6y, 5^, and TOq that modify the masses of the flavor- 
neutral physical-basis pions by a factor of 1/4. Explicit factors of these parameters 
remain unchanged, and this arises from the fact that the implicit factors have a 
one-to-one correspondence with the internal quark loops. The former case, which 
corresponds to the connected terms, is also rather straightforward. Identifying the 
quark loops leads to the conclusion that they must each be multiplied by a factor 
of 1/4. This factor will be directly seen in our results unlike the correction to the 
disconnected pieces. More detail on quark flow analysis for the specific calculations 
to follow can be found in Ref. [551 1311 1331 

At this point, we can calculate directly from the Lagrangian, Eq. p8)) . quantities 
such as one-loop masses and decay constants, or pion scattering amplitudes. Of 
these, we will discuss the masses and decay constants. We can also calculate weak 
matrix elements, such as Bk, although this calculation is much more involved, so 
we will refer the reader to Ref. [34|for more details. 



4.1. Pion masses 

To calculate the pion mass to one loop in SXPT, we need to expand out the La- 
grangian to quartic order in $ and use the four-point vertices that arise to calculate 
the diagrams shown in Fig. [21 These are the contributions to the self-energy of 
the pion field, and the first diagram includes those which are connected, meaning 
they come from the l/(q^ -I- m^) parts of propagators (both fiavor-charged and 
fiavor- neutral) . The second diagram with the cross in the propagator includes those 
diagrams with only the flavor- neutral axial, singlet, and vector taste disconnected 
propagator, T). 



I 



I 
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(a) 



(b) 



Fig. 2. The two one-loop diagrams contributing to the self-energy of the pion, and thus to its 
mass, (a) includes the connected diagrams and (b) includes the disconnected diagrams. 



We will also treat our calculation as a partially-quenched theory, where some 
quarks are treated as valence quarks while others are sea quarks. The valence quarks 
will not be allowed to propagate in loops. By using the same quark flow analysis to 
correct for the rooting trick, we can take these diagrams out by hand and not include 
them in our result. This is done in detail in Ref. |2S]for the pion masses, so we will 
not discuss this further here. The more formally correct method of partial que nching 
is to upgrade our symmetry group from SU{4:N) to 5L/(4iVsca + 4iVvai|4iVvai),'^^'^ 
a graded group which can explicitly remove valence quark loop diagrams. For the 
cases discussed throughout this review, the quark loop method is simpler to employ, 
and leads to identical results. 

To show only results that are specifically relevant to lattice simulations, we will 
pick TV = 5, where there are two valence quarks with masses rrix and my, and three 
sea quarks with masses mi, mi, and m^. The up and down sea quarks are set to be 
degenerate, and the strange quark mass will be allowed to have a different mass. 
Picking out the "real" pion amounts to setting m^ = my = mi while the "real" 
kaon requires us to set m^ = mi and my = mg. Also, we will only calculate the one- 
loop Goldstone pion, or the pseudoscalar-taste pion, as this is the pion whose mass 
vanishes in the chiral limit at finite lattice spacing (recall for this case, Ap = 0). 

The details of the calculation can be found in Ref. [25(, and the result, including 
the tree-level analytic terms discussed above is 



l-loop\ ' 



fx{mx + my) 



167r2/2 



^ E (-^3; Mi) l{mlj) -f 1^(2^8 - L5)(m, -h m,) 



(27) 



with Li the Gasser-Leutwyler low-energy constants,'^ C is an unknown constant 
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which comes from terms in £4 that involve the lattice spacingpSl and 

i{m') = rn' In (^^ j , Rf'"^ {{m}■,{^^}) ^ ^ , (28) 

J|(to- -TOj) 

i=l 

where A is the chiral cutoff scale. The sets of masses given in the residues R are 
those appearing in the numerators and denominators of the disconnected parts of 
the flavor-neutral propagators. More explicitly: 

In Eq. ([27]) , the sums over j run over the masses in the first argument of the residue 
functions R. 

We will discuss the results from using this expression to fit lattice data after we 
calculate the decay constant. 



4.2. Pion decay constants 

The diagrams needed to calculate the pion decay constant, again in the partially- 
quenched case, are shown in Fig. [3l where the notation is the same as before. The 
addition here is the solid box, which represents an insertion of the axial current 
corresponding to the pseudoscalar taste pion, given at the chiral level by 



■Xy,P rp 

Ja<5 = Tr 



(29) 



with 7-'^'' a projector which extracts the 4x4 block matrix which corresponds to 
the charged pion.12^ We can extract the decay constant from the matrix element 



■xy,P 



(30) 



with P5 the pion state with quark content xy. With this normalization, w 
131 MeV. Note the correction coming from the wavefunction rcnormalization 
(and thus from the self-energy graphs above), are proportional to the vertex 
corrections .^221 

The result for the decay constant is much more complicated than that of the 
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(a) (b) (c) (d) 



Fig. 3. The four one-loop diagrams contributing to the decay constant of the pion, using the 
same notation as Fig. [2] The box corresponds to an insertion of the axial current, (a) and (b) 
are the wavefunction renormalization diagrams (connected and disconnected, respectively) and (c) 
and (d) are the vertex corrections 



/"l-loop 
J xy. h 



= 1 



V i 

+dly E 4'}" (-^5; Mi) £{ml) + 2 E {M^ Mi) l{ml)\ 
j j ' 

+ \a'5'v f C E i^^' ^2) eimlv) 

V j 

+9L E i^r; MX) l{mly) - 2 E <y'' {M[, Mi) €(™,V 
j 



+ 



{2me + rus) L4 + {m^ + ruy) L5 + a^F , (31) 



where we have used the shorthand notation = d/drn^^ j, and similarly for the 
other tastes and flavors. We have the additional mass sets 

Ml = {mxx,umr^^t} , Ml = {niyy^t.rnrj.t} , 
Ml = {mxx,t,mrj.t,mr,',t} , Mj = {myy^t,mrj.t,rn^' .t} ■ 

Again, the sums over j run over the mesons in the sets which are the first argument 
of the corresponding i?'s. The sum over t is over the 16 tastes and the sum over q 
is over the sea quarks u, d, and s. 



4.3. Results from pion fits 

The MILC collaboration has fit to Eqs. (P?]) and (|3ip with multiple lattice spacings 
(a « 0.125 fm, 0.09 fm, and 0.06 fm). As one can see from Refs. [HI 133 the fits are 
quite good, with a confidence level of 0.99. There are two subsets of fits, described 
in Ref. [SI giving either 122 or 978 data points, and so the fits are constrained 
considerably. The most current results (still marked as preliminary), after taking 
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the physical quark mass hmit and the continuum limit are, for the decay constants 
= 128.6 ± 0.4 ± 3.0 MeV , = 155.3 ± 0.4 ± 3.1 MeV , 



where the first errors are statistical and the second systematic, and for the quark 
masses (adding all the errors in quadrature) 



One can also extract the Gasser-Leutwyler coefficients, the up and down quark 
masses, and certain CKM matrix elements, but we refer the reader to Ref . [6l [37l for 
those quantities. 

The key point from these results is that with Eqs. ((27)) and (|3T|) . one can ob- 
tain rather precise results for certain quantities of interest with staggered fermions. 
Without these expressions, one would have no information about what the true chi- 
ral behavior of the pion mass and decay constant is near the continuum limit. In fact, 
for comparison, the confidence level of fits when using standard chiral expression92l 
[i.e., the continuum limit of Eqs. ([27]) and ([3T|) ] is ^ lO^^'"^". Using these expressions 
is crucial for a result with small, credible, errors. 

5. Heavy-lights 

Having had success describing the simplest of quantities, we now turn our attention 
to more difficult quantities which are equally important for phenomenology. That is, 
including heavy-light mesons in our SXPT. With this, we can determine the chiral 
behavior of heavy-light decay constants, semileptonic form factors, as well as m ore 
complicated quantities such has B Div or B D*£v form factors . ^^^^^ ^ ^^ ^ 

The procedure for including heavy-light mesons is rather straightforward. We 
will sketch a few steps and refer the reader to Ref. [iPl HTl [551 for more details, both 
on the continuum and staggered version of this procedure. After doing this, we will 
calculate the heavy-light decay constant for a B meson. 

The first step is incorporate heavy quark effective theory (HQET) into SXPT. 
We take the masses of the heavy quarks to be large (compared to Kqcd) and 
perform a l/mq expansion. For this review, we will only keep leading order terms 
in this expansion [0(1) terms]. A heavy-light field, if^, has two indices, one we'll 
denote as a superscript Q to refer to the heavy quark flavor and the other will be 
a to denote the light quark flavor/taste. We also assume that the heavy quark is 
discretized in such a way that the doublers have large enough masses so they don't 
affect the dynamics (see Ref. [33|for a more detailed discussion of this restriction), 
and so we treat it as a continuum quark. 

For m heavy quarks, we have a combined spin-flavor SU (2m) symmetry in ad- 
dition to the chiral SU{AN)l x S'/7(4iV)fl symmetry. In matrix notation, we have 



/a7/.- 1.208(2) , 



mf s = 90(6) MeV , mf^ = 3 
ms/mi = 27.2(4) . 



3.3(2) MeV , 



(32) 
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where S G SU{2m), and U e SU{4:N). We pick a special transformation for the 
right-hand side, so as to keep iJ's transformation under parity simple. This matrix 
U is spacetime dependent and is defined with cr^ = E, where a transforms under 
the chiral symmetry group as 

a LcrU^ = VaR . (33) 

When we enumerate all of the possible operators that is part of the heavy- light 
SXPT action, we find there are over 200 terms in addition to the continuum action! 
All of these additional terms, however, contribute in the same way only to analytic 
terms when working to one-loop order in any given quantity. As these terms are all 
proportional to a^, however, they are unphysical quantities; thus it is not necessary 
to determine precise values for them all. For the following calculation, we will just 
write down the relevant action required for the one-loop terms which involve the 
heavy- lights: 

Chi = Tr (hHv -D^+g^ Tr (HHj^'j^A^) , (34) 

where the trace is a combined trace over the heavy quark indices, the light quark 
flavor-taste indices and the Dirac indices, and 

^ - d^Ha + iHbWf , (35) 
Y^='-[cT^d^,a + ad^a^] , (36) 

A^ = ^[c7^d^,a^ad^a^] . (37) 

One can verify that Eq. ([M]) is invariant under the full HQET-SXPT symmetry 
groupB There are terms which arise due to the finite light quark masses and the 
lattice spacings, but these will only contribute to analytic terms, and so we will not 
include them here. 

To calculate Jb, we need the left-handed current. To leading order, this takes 
the form 

r^ = !^tvDh^{l-75)H)a^X^^^ , (38) 

where A^'^' is a constant vector which fixes the taste-fiavor: (A'"^-*)^ = Sac and tr^i is 
a trace over only Dirac indices. We have not included an index on the heavy- light 
field for the heavy quark content, although below we will refer to a B meson, a 
heavy-light made from a bottom quark and a light quark, to stress the fact that 
we are working to leading order in the heavy quark mass. The decay constant is 
defined through the matrix element 

(0|j^'=|B,(w)) = ifB^niB^v^Sac , (39) 

®Note that Eq. I I34II is in Minkowski space, because it is more consistent with other (continuum) 
work with heavy-lights. By changing the signs of the mass term and V in Eq. Ijl8|l . one can obtain 
the SXPT action in Minkowski space as well. 
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(a) 



(b) 



(c) 



(d) 



Fig. 4. The four one-loop diagrams contributing to the B decay constant, (a) and (b) are the 
wavefunction renormahzation diagrams (connected and disconnected, respectively) and (c) and 
(d) arc the vertex corrections. 



where the normahzation of the B state is given in Ref. 1331 At leading order we 
have f]^ = nj ^niBa ■ The non- vanishing diagrams are shown in Fig. 31 where the 
notation is the same as in the pion case, except we now have heavy-Ught mesons, 
denoted by the double lines. We need to account for the factors of 1 /4 to correct for 
the fourth-root trick, and it turns out through a careful quark-flow analysiJ^Sl to 
arise in the same way as for pion quantities. The results below have already taken 
this into account. 

We have for the partially quenched decay constant for 2+1 flavors 



/-Bx 
fLO 



= 1 



-I- Cs{2mi + TUs) + c^m^ + CqC^ 



(40) 



where the sums over i, g, and j are the same as in the pion case, and the mass 
sets are defined above. The analytic terms are combinations of numerous unknown 
parameters from the Lagrangian as well as higher-order terms in the current. Finally, 
due to a shift symmetry in the HQET-SXPT Lagrangian, this result is independent 
of the taste of the light quark in the B meson.l^ 

Although we have this calculation in terms of the B meson, in principle it is also 
valid to leading order for the D meson, although some errors may be larger due to 
the smaller mass of the charm quark. Nevertheless, the D and Ds decay constants 
have been measured by the Fermilab Lattice group!^ 

/c+ = 201(3)(17) MeV , fo^ - 249(3)(16) MeV , (41) 

which agrees with the experimental measurement that was posted shortly 
afterwards.S^ Again, from these results, one can extract CKM matrix elements.^ 

Using this HQET-SXPT theory, we can also calculate the semileptonic form 
factors for B — > -kHv decay, B D*£i', or many other quantities one might be 
interested in calculating on the lattice. We do not discuss these results in detail, 
but instead refer the reader to Ref. [351 1331 ISSl for details. 
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6. Conclusions 

We have discussed many of the interesting applications of SXPT to extract physical 
results from lattice data. There are many other quantities we have not discussed, but 
would like to touch upon here. Van de Water & Sharpe have analyzed the neutral 
kaon mixing parameter, Bk, using SXPT.'S^l This is rather extensive and requires 
a large number of additional parameters to control all of the errors that arise in a 
staggered calculation of Bk (operator mixing, perturbative errors, etc.). Perhaps 
this is an indication of the practical limitations of staggered quarks, as when the 
quantities become more complicated (as is the case with weak matrix elements), 
the difficulty increases dramatically. An interesting solution for Bk and other weak 
matrix elements could be a "m ixed-actio n" approach, using different valence quarks 
on top of staggered sea quarks .SSEUli! 

A slightly different application is to that of a calculation of the muon g — 
2| 46 | 47 | 48|49 | ^ calculation of the muon g — 2 on the lattice is actually a calculation 
of the photon vacuum polarization, and using SXPT combined with electromagnetic 
interactions for this can aid in extracting the low-energy behavior of this function. 
For this, an extension of SXPT to include light resonances (mainly the vectors) is 
needed, as they dominate the low-energy behavior. Although this has been used suc- 
cessfully to describe the lattice results, the errors are not yet competitive enough 
with other techniques to calculate the muon g — 2 theoretically (see Ref. |50] for 
example). 

In conclusion, at least for the simpler quantities, precise staggered simulations 
combined with SXPT expressions allow for accurate determinations of physical 
quantities that can be compared against experiment reliably. This is an indica- 
tion that lattice QCD can be used in the forthcoming years as a spectacular testing 
ground for more precise determinations of the effects of the strong interactions. 
Whether or not more complicated quantities can be adequately determined using 
staggered quarks is an open question, but for the present, the use of SXPT has made 
possible these precise extrapolations. 
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